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Abstract 

The two-orbital J\-Ji model that describes a square lattice of frustrated spin-1 iron atoms is 
analyzed within the linear spin-wave approximation and by exact diagonalization over a 4 x 4 
cluster. A quantum critical point (QCP) is identified that separates hidden magnetic order at 
weak Hund's rule coupling from a commensurate spin density wave (cSDW) at strong Hund's rule 
coupling. Although the moment for cSDW order is small at the QCP, the critical linear spin-wave 
spectrum shows strong low-energy excitations centered at the wavenumbers that correspond to 
cSDW order. These disperse anisotropically. A fit to the magnetic excitation spectrum of ferro- 
pnictide materials obtained recently by inelastic neutron scattering measurements notably accounts 
for the absence of softening at the wavenumber that corresponds to Neel order. 
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I. INTRODUCTION 



The discovery of ferro-pnictide superconductors has set a new and unforseen direction in 
the search for high-temperature superconductivity^ A common structural feature of these 
materials are square lattices of iron atoms that are stacked up on top of each other.- Density 
functional theory (DFT) calculations find that the iron 3d orbitals contribute most of the 
electronic spectral weight near the Fermi surface.- 1 ^ This implies that the iron 3d electrons 
surprisingly play an active role in the superconductivity. In contrast, DFT calculations also 
predict that the iron moments order antiferromagnetically along the a and along the c axes 
of the iron lattice at low temperature in the undoped parent compounds, 5 with an ordered 
moment of approximately 2 Bohr magnetons {hb)-~ Elastic neutron diffraction studies have 
confirmed the prediction of long-range commensurate spin-density wave (cSDW) order, but 
they find a much lower ordered moment that can be a small fraction of the Bohr magneton 
(0.36 /is) in some cases. 2 It was suggested early on that Heisenberg exchange between both 
nearest-neighbor ( Ji) and next-nearest neighbor ( J 2 ) spins on a square lattice of iron atoms 
could account for the low ordered moment that is observed because of magnetic frustration.- 
Recent calculations based on the linear spin-wave approximation 8 indicate that this effect 
requires fine tuning however,- and that it occurs only very close to the point of maximum 
frustration, J2 = 0.5 J\. A possibly more viable explanation is that the low ordered cSDW 
moment commonly observed in ferro-pnictide parent compounds is due to proximity to 
hidden magnetic order among particular iron 3d orbitals such that the magnetic moments 
cancel out per iron atom. The author and Rezayi have shown that such violations of Hund's 
rule can occur in multi-orbital Heisenberg models with off-diagonal frustration.— 

Spin-wave excitations of the cSDW groundstate have also been detected in parent 
compounds to ferro-pnictide superconductors at low temperature by inelastic neutron 
scattering.—^ 2 - Contrary to the prediction of the conventional J1-J2 model over the square 
lattice however,— the measured spin-wave spectrum fails to soften at the wave number as- 
sociated with Neel order, (7r/a,7r/a). Here a denotes the square lattice constant. It shows a 
local maximum there instead.— Fits to the measured spin-wave dispersion invoke unphysi- 
cally large anisotropy in the Heisenberg spin exchange across nearest-neighbor links on the 
square lattice of iron atoms in order to account for the dispersion at the Brillouin zone 
boundary. More recently, low-energy magnetic excitations centered at the cSDW wave num- 
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bers have been detected in a ferro-pnictide superconductor.^ These disperse up in energy in 
an anisotropic fashion that is consistent with the presence of some degree of cSDW order. 

In this paper, we continue to explore the hidden magnetic order that is predicted to exist 
in multi-orbital Heisenberg models over the square lattice when the diagonal and the off- 
diagonal exchange coupling constants are sufficiently different.— Two spin-1/2 orbitals per 
iron atom are assumed, which implies spin-1 iron moments if Hund's rule is obeyed. A quan- 
tum critical point^ is identified within the linear spin-wave approximation that separates 
hidden ferromagnetic or hidden Neel order at weak Hund's rule coupling from a cSDW at 
strong Hund's rule coupling. Numerical exact diagonalizations of the two-orbital Heisenberg 
model over a 4 x 4 lattice indicate that the moment associated with cSDW order decays 
strongly as it enters hidden magnetic order through the putative quantum critical point. 
The linear spin-wave approximation, on the other hand, finds observable low-energy spin- 
wave excitations at the quantum critical point that disperse linearly from the cSDW wave 
numbers in an anisotropic fashion. A spin gap begins to grow there upon entering hidden 
magnetic order. Observable low-energy spin-wave excitations due to hidden order that are 
centered at zero momentum are also predicted by the linear spin-wave approximation, but 
these have very small spectral weight. Fits of the spin- wave spectrum observed by inelastic 
neutron scattering in ferro-pnictide material a ^ ^ to the critical linear spin-wave spectrum 
are achieved. They notably account for the dispersion of the low-energy spinwaves centered 
at the cSDW wave numbers and for the local maximum that exists at the wavenumber 
associated with Neel order. A critical Hund's rule coupling of order 100 meV is extracted 
from these fits. Below we discuss useful theoretical aspects of the Heisenberg model for 
ferro-pnictide materials 7 before we go on to analyze it. 

II. MULTI-ORBITAL HEISENBERG MODEL 

Recent low-temperature measurements of the optical conductivity in the cSDW phase 
of parent compounds to ferro-pnictide high-T c superconductors find a strong suppression of 
the integrated kinetic energy for the conduction electrons compared to that predicted by 
electronic band-structure calculations.— A low carrier density in the cSD W phase of ferro- 
pnictide materials is consistent with quantum-oscillation experiment a l ^ and with angle- 
resolved photoemission studies^ that find evidence for tiny 2D Fermi-surface "spots" , each 
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on the order of 1% of the Brillouin zone. Also, the dependence on temperature shown by 
the electrical resistance of such parent compounds indicates proximity to a metal-insulator 
transition.— Undoped ferro-pnictide materials in the cSDW phase then may very well lie 
near a Mott transition, where conduction electrons experience moderately strong repulsive 
interactions. By continuity with the limit of strong inter-electron repulsion, we therefore 
believe that a local-moment description of magnetism in parent compounds to iron-based 
high-T c superconductors is a justifiable starting point at low temperature. The absence of a 
Stoner continuum of incoherent magnetic excitations reported in recent experimental studies 
of magnetic excitations in the cSDW phase of a ferro-pnictide material is consistent with this 
view. 14 Consider then the following multi-orbital spin- 1/2 Hamiltonian that contains near- 
neighbor Heisenberg exchange among local iron moments within isolated layers in addition 
to Hund's-rule coupling: 7 

^ = ^oEfEs,(«)| 2 + EE^' /3 s l («)-s,(/3)+ e E^m-s^). (i) 

i 1 a J (hj) a£ «tj» a,/3 

Above, Si (a) is the spin operator that acts on the spin- 1/2 state of orbital a in the iron atom 
at site i. The latter runs over the square lattice of iron atoms that make up an isolated layer. 
The application of Hund's rule is controlled by a local ferromagnetic Heisenberg exchange 
constant J < 0, while nearest neighbor and next-nearest neighbor Heisenberg exchange 
across the links and is controlled by the tensor exchange constants Jf and 

J 2 i respectively. The strength of the crystal field at each iron atom determines the number 
of orbitals per iron atom above. It can be as low as two for strong crystal fields compared 
to the Hund's-rule coupling, in which case spin-1/2 moments exist on the iron 3d xz and 
3d yz orbitals.- 1 ^ Four 3d orbitals per iron atom carry spin-1/2 moments in the case of weak 
crystal fields compared to the Hund's-rule coupling, on the other hand. 3,7 Last, the above 
model Hamiltonian reduces to the conventional J\-Ji model over the square lattice when 
Hund's rule is enforced,-"^ —Jo oo, where the spin at each site is given by half the 
number of orbitals, S — 1 or S — 2. 

We now illustrate a useful transformation of the model Hamiltonians (CQ) over the square 
lattice that relates eigenstates and eigenvalues with different nearest-neighbor Heisenberg 
exchange coupling constants. Let p denote a permutation of the orbitals. Consider then the 
partial "reflection" Pa{b)'- 

orbital a — > p(a) on the A(B) sublattice (2) 



of the Neel state. If \I/ is an eigenstate of the original model Hamiltonian ([I]) with energy 
E, the transformed state = Pa(b)^ is clearly then an eigenstate of the transformed 
Hamiltonian H' = Pa(b)HP^ b ^ with energy E. Importantly, the Heisenberg exchange 
coupling constants transform under ([2]) as 

when orbital a lies on the A{B) sublattice of the Neel state. Both J%'^ on the B(A) sublattice 
and the Hund's rule coupling J remain fixed. Yet how does the spectrum of states transform 
under ([2])? Let T x and T y denote translations along the x axis and along the y axis by one 
lattice constant a. Because T x T y moves each sublattice to itself, it commutes with Pa and 
with Pb- If \P has a crystal momentum hh, we therefore have that T x T y ^' = e l ( k ^+ k y) a v]/'. 
The crystal momentum then transforms under (TS]) as 

k — > k or k + (ir/a, ir/a). (4) 

Also, the spin of the new state ^' remains unchanged because Pa and Pb both commute 
with Si = J2 a Si(a). 

It is also useful to consider the special case where all nearest-neighbor and next-nearest- 
neighbor exchange coupling constants are equal, respectively: J" = J\ and J^ ,/3 = J<i- The 
model Hamiltonian ([1]) then reduces to H = | J J2i Si • Sj + J\ Yluj) Si • Sj + J2 Ylaij)) Si • Sj. 
Because Si + Sj commutes with Sj • Si, the latter then commutes with the Hamiltonian. 
The total spin at a given site i is hence a good quantum number. This means that the 
groundstate obeys Hund's rule in the classical limit because states with maximum total spin 
at a given site minimize both the Hund's-rule energy ( J < 0) and the Heisenberg exchange 
energies in such a case. A violation of Hund's rule will therefore require some variation in 
the Heisenberg exchange coupling constants among the different iron orbitals. 

III. SPIN-1 IRON ATOMS 

The iron 3d orbitals in ferro-pnictide materials experience a crystal field due to the neig- 
boring pnictide atoms that lie above/below the center of each unit square of iron atoms. Both 
DFT calculations and basic considerations indicate that the crystal field splits these levels so 
that the degenerate 3d xz and 3d yz orbitals lie just below the highest-energy 3c? orbitaL 4,22 The 
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Pauli principle then implies that spin s = h/2 moments exists on both the 3d xz and the 3d yz 
orbitals in the limit of weak Hund's rule coupling. Henceforth, we shall assume that only two 
spin-1/2 orbitals, a and b, exist per iron atom in parent compounds to ferro-pnictide high-T c 
superconductors. This is consistent with the ordered moment of 2/i# obtained by DFT 
calculations for antiferromagnetic groundstates of ferro-pnictide parent compounds.-^ Let 
us next fix notation by assuming unique diagonal and unique off-diagonal Heisenberg ex- 
change coupling constants over the nearest neighbor and the next-nearest neighbor links of 
the square lattice of iron atoms: 

ja,a t[| jb,b 1 ja,b t_|_ jb,a 

J l(2) — J l(2) — J l(2) ariCl J l(2) — J l(2) — J l(2)- l D J 

The author and Rezayi have recently shown that strong enough off-diagonal frustration, 
J 2 > 0, can lead to either hidden ferromagnetic order or to hidden Neel order among spin- 
1/2 moments that are governed by the J0-J1-J2 model Hamiltonian= ([I]). Figure [T] displays 
the two types of hidden order, both of which show no net magnetic moment. Proximity of the 
collinear antiferromagnetic state, which exists at J2 > \Ji\/2 when Hund's rule is enforced, 
to such hidden order can therefore account for the low ordered moment that is observed 
by elastic neutron diffraction in the cSDW state that exists at low temperature in undoped 
ferro-pnictide materials.^ Balancing the energy of such hidden-order groundstates in the 
classical limit, s — > 00, against that of conventional magnetic order leads to critical values 
for the Hund's rule coupling that are listed in Table [H True magnetic order becomes unstable 
to hidden magnetic order at weaker Hund's-rule coupling, — J < — Joe- A true magnetically 
ordered groundstate therefore first becomes unstable to the hidden-order groundstate that 
has the larger of the two critical values for the Hund's rule coupling strength, — Jq c . By 
Table QJ for example, the true cSDW state becomes unstable to hidden ferromagnetic order 
at weak Hund's rule coupling only if the off-diagonal nearest-neighbor Heisenberg exchange 
is larger than the corresponding diagonal exchange: J-j 1 > jj . 

Linear Spin- Wave Theory. Study of the classical limit, s — > 00, therefore indicates that 
extremely low ordered moments are possible at weak enough Hund's rule coupling, J < 0, 
when off-diagonal frustration exists: J 2 > 0. The nature of spin- wave excitations about the 
hidden magnetic order that is shown in Fig. [Dean also be obtained from the classical limit. 
In such case, each spin-1/2 moment obeys the following dynamical equation for precession: 

x Si(a). (6) 



la 



OH 
dS t (a) 



a 
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HIDDEN FERROMAGNETIC ORDER 



HIDDEN NEEL ORDER 




FIG. 1: Shown are two groundstates of a square lattice of spin-1 iron atoms that violate Hund's 
rule, but that exhibit hidden magnetic order. [See Eq. ([!]) and Table HI] 





True cSDW 


True Neel 


True Ferromagnet 


Hidden Ferromagnet 
Hidden Neel 


-Joc = 2(J r - J}) - 4jJ 
- Joc = 2(j| - J r ) - 4j| 


- Jqc = ^2 ~ 4J 1 

- Joc = ±J-t - 4Jr 


- Joc = Uj- + 4J^ 

- J 0c = \J^ + 4J| 



TABLE I: Listed are the phase boundaries that separate true magnetic order from hidden magnetic 
order. The critical Hund's rule coupling is obtained by balancing the classical energies of the model 
Hamiltonian JT]). 

Here we have applied a magnetic field to the system, hj(a), that is in principle sensitive to 
the orbital degree of freedom, a = a or b. If the external magnetic field is a plane-wave with 
wavenumbers k and frequency u, then its transverse component with respect to the spin 
axis for hidden magnetic order generates a corresponding linear response in the transverse 
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components of the electronic spin. It has the form 



S ± (±, k, u) = x±(±, k, co) h ± (±, k, u) (7) 

and 

Sj_(=f, k, lo) = ix±(K u) ny x h ± (±, k, to) (8) 

in the case of the hidden ferromagnetic state. Here we define true (+) and hidden (— ) 
magnetic moments and magnetic fields by S(±) = S(a)±S(6) andbyh(±) = [h(a)±h(6)]/2, 
while n|| denotes the unit vector parallel to the magnetization axis for hidden order. Recall 
that the hidden ferromagnet is actually an antiferromagnet over the two sublattices set by 
the two orbitals a and b (see Fig. [T]and ref.— ). The true magnetic susceptibility, k, uj), 

the hidden magnetic susceptibility, x_i_( — > k, u), and the conjugate susceptibility x_i_(k, oj) 
can then be obtained by solving the linear precession equation (jH]) that results in terms of 
the four transverse spin components per site, Sj_(a) and Sj_(6). This yields 

X±(±, k, w) = s J^TyT2^ + i^-^T 1 - [« - (^-) V T ') (9) 

and 

X±(K") = /2 ([u + (fi+fi-) 1 / 2 ]- 1 - [w - (n+O-) 1 / 2 ]- 1 ), (10) 

with frequencies given by 

= s(Jt-4) "£(2sm~k n a) 2 + s(j£-4) £ (2sin^ n a) 2 

n=x,y ^ n=+,— 

Q+ = 2sJq + sJi ^ (2 cos -k n a) 2 — sjf V] (2sin-/c n a) 2 

n=x,y ^ n=x,y ^ 

+sJ^ 51 (2cos^A;„a) 2 - sjj ^ (2 sin ^„a) 2 

n=+,- ^ n=+,- ^ 

in the case of the hidden ferromagnetic state. Above, k± — k x ± a denotes the square 
lattice constant, and s is the electron spin h/2. Next observe that the hidden Neel state is 
related to the hidden ferromagnetic state by the transformation (j2J) mentioned in the previous 
section. (See Fig. [U) By ([3]), the frequencies 0± in that case are therefore obtained from 
the ones listed above for the hidden ferromagnet after making the exchange of the nearest- 
neighbor Heisenberg coupling constants, jj <H- Jj 1 . Further, the wavenumbers of the hidden 
spin S(— ) and of the hidden magnetic field h(— ) that appear in the linear-response equations 
((7|) and (jSJ) must be replaced by the wave numbers on the antiferromagnetic sublattice, k — >• 
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k— (ir /a, ir /a). The wavenumbers of the true spin and of the true magnetic field that appear 
there remain unchanged, however. The end result is that the true magnetic susceptibility 
X_i_(+, k, u) in the hidden Neel state is given by that in the hidden ferromagnetic phase with 
the exchange j\ -H- J-j 1 alone, whereas the hidden magnetic susceptibility x±(~ )k, u)) there 
is given by that in the hidden ferromagnetic phase with the previous exchange along with 
the replacement k — > k + (it /a, ir/a) in the wave numbers. [Cf. Eq. (HJ.] 

Equation (jUJ) and the discussion above yield true static magnetic susceptibilities x_l(0) = 
(J + 47-j 1 + ^J^y 1 for the hidden ferromagnetic state and x±(fy = (Jo + 4 4' + ^Jt)~ l 
for the hidden Neel state. It also yields static susceptibilities for hidden magnetic order 
given by ,k, 0) = 2s/Q_, which correctly diverge as k approaches the wavenumber 
that corresponds to either long-range ferromagnetic or Neel order. Finally, the poles of 
the dynamical susceptibility for hidden magnetic order fl9]), x±(~ > k, uj), represent two spin- 
wave excitations of momentum frk that correspond to the two possible orientations of the 
transverse component of the hidden moment, Sj_(— ,k, ou). They both lie at an energy 
^(fi+fi-) 1 / 2 above the groundstate energy. These spin-wave excitations are also reflected 
by the poles in the dynamical susceptibility for true magnetic order (jHJ), k, u). As 

mentioned above, however, the wave number k shifts by the characteristic amount (ir/a, ir/a) 
in the case of hidden Neel order. Last, both spinwaves have spectral weights 

A(±,k) = 7cs{Q T /{l±) 1/2 (11) 

in the true (+) and the hidden (— ) channels that are read off from the imaginary part of 
the dynamical susceptibility: 

Imxi(±,k,u) = A(±,k)(8[u- (Q + tt_) 1/2 } -5[u + (Q + Q_) 1/2 }). (12) 

The peaks in the dynamical structure function above signal spin-wave excitations that are 
observable by inelastic neutron scattering only in the true magnetic (+) channel. Recalling 
the discussion following Eq. (jUJ), the true dynamical susceptibility X-i-(+> k, oj) is common to 
a hidden ferromagnet and to a hidden Neel state that are related by the transformation 
Inelastic neutron scattering is therefore unable to discriminate between these two possible 
types of hidden magnetic order! 

Figure [2] depicts the linear spin-wave spectrum for a cSDW with intervening hidden 
ferromagnetic (Neel) order at two extreme values for the Hund's rule coupling: absent and 
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SPIN-WAVE EXCITATION ENERGY IN ABSENCE OF HUND'S RULE COUPLING 




SPIN-WAVE EXCITATION ENERGY AT TRANSITION INTO HIDDEN MAGNETIC ORDER 




FIG. 2: Displayed is the linear spin-wave spectrum of the Hamiltonian ([I]) for the following set 
of Heisenberg exchange coupling constants: jf = 0, J^~^ > 0, and jj = 0.3 J^~^ = J^~. The 
Hund's rule coupling at the quantum critical point that separates hidden ferromagnetic (Neel) 
order from a cSDW is set by the value listed in Table [H Spin-wave energies are given in units 
of fi 2 J± and the wavenumber k' is associated with the antiferromagnetic sublattice for hidden 
order. (See Fig. [TJ) 

critical. Here, we have chosen Heisenberg exchange coupling constants J^ ±S) = 0, > 0, 
and jj = 0.3 J{ L ^"' ) = J^~, which lies near maximum frustration when Hund's rule is obeyed. 
The top graph displays how the spectrum collapses linearly to the groundstate energy at 
zero momentum in the absence of Hund's rule coupling. In general, the spin-wave velocity 
is given by 

v = 2sa{[Ji - j\ + 2(4 - jj)] • [i J + 2Ji + 2Ji}) 1 ' 2 (13) 

in the hidden ferromagnet, and by the same expression after the exchange (j2J) j\ J\ in 
the hidden Neel state. It correctly collapses to zero at the phase boundary with the true 
ferromagnetic state by Table HI at which point the spin- wave excitation frequency disperses 
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(0,0) (ji,0) (it.Jt) (0,0) 



MOMENTUM (1/a) 

FIG. 3: The spectral weight of observable (+) spinwaves at the quantum critical point that sepa- 
rates hidden ferromagnetic (Neel) order from a cSDW is shown in units of H. The following set of 
Heisenberg exchange coupling constants are used: = 0, J{^"^ > 0, and jj = 0.3 J-j 1 = J^~. 

The critical Hund's rule exchange coupling is set to the value listed in Table HI Also shown are 
contours over the Brillouin zone for such spinwaves at a low fixed excitation energy. 

quadratically with wavenumber. Notice, however, the mild depression in the top spectrum 
that is centered at the cSDW wave number \s. c sdw = (tt/o, 0). The bottom graph in Fig. [2] 
shows how it softens completely to zero energy at the critical value for Hund's rule coupling, 
— Jq c = 0.8 Ji^\ in which case the hidden ferromagnetic (Neel) state transits into the cSDW 
state by Table HI These special zero-energy modes result from the collapse of the frequency 
Q + (k c SDw) to zero there. Evaluating (fi+fi.) 1 / 2 at \l c sdw from the expressions for the 
frequencies that are listed below Eq. f flOl) yields the spin gap there: 

A cSDW = 2s[(4J 2 x - J 0c )(J - Joc)] 1/2 - (14) 
The former also implies a divergent spectral weight (|lip for such spinwaves at the quantum 



11 



phase transition: 

A(+, KSDW) = 7TS[(U£ ~ J c)/(Jo ~ J c)} 1/2 - (15) 

Because inelastic neutron scattering is proportional to lmx±{+, k, u), low-energy spin-waves 
with momentum near Kk c sDW wm dominate such measurements near the transition between 
the cSDW state and hidden magnetic order (see Fig. [3]). Study of the expressions for the 
frequencies Q± that are listed below Eq. (11 01) yields that the dominant spin-wave excitations 
centered at the cSDW wave numbers \t c SDW disperse up from the groundstate energy as 
[vf (k[ — ir /a) 2 + v 2 k 2 + A' 2 SDW ] 1 ^ 2 in general. (See Fig. [2j) Here k\ and k t are the components 
of k that are respectively parallel and perpendicular to \h c sdw- Near criticality, A cSDW — > 0, 
the cSDW order is short range, with correlations lengths ^sdw = v i/^cSdw an d £,csdw = 
Vt/ A c sdw m the longitudinal and transverse directions. At criticality, A c sdw = 0, the 
longitudinal velocity v\ coincides with the spin-wave velocity u for hidden magnetic order 
f|T3|) . while the transverse velocity v t is equal to the previous divided by the anisotropy 
parameter 

vi 



2(4 + jj-) + 4 + j x ±iV2 



(16) 



2(4 + Jir) ~ 4 ~ Ji 

The latter coincides with the anisotropy of low-energy spinwaves that is predicted by the 
conventional J\-Ji model, in which case Ji( 2 ) = [4(2) + *A(2)]/2 by Hund's rule. 

Exact Diagonalization. The above semi-classical analysis indicates that large enough off- 
diagonal frustration can induce a quantum phase transition into hidden magnetic order that 
is unfrustrated, but that violates Hund's rule. We shall now confirm this by obtaining the 
exact low-energy spectrum of states for the J0-J1-J2 model (CQ) with two spin- 1/2 moments 
per site over a 4 by 4 square lattice with periodic boundary conditions. The model Hamilto- 
nian operator (JTJ acts on a Hilbert space that is restricted to have 16 up spins and 16 down 
spins. Next, translational invariance is exploited in order to bring H into block-diagonal 
form. Each block of the Hamiltonian is labeled by the allowed momentum quantum num- 
bers: (k x a, k y d) = (0, 0), (tt, 0), (it, it), (vr/2, 0), (tt/2, ir/2) and (tt, 7r/2), plus their symmetric 
counterparts. Spin-flip symmetry is also exploited to further block-diagonalize the Hamil- 
tonian at such momenta into two blocks that are respectively even and odd under it. The 
combination of translation and spin-flip symmetries reduces the dimension of each block to a 
little under 19,000,000 states. Next, we construct the Hamiltonian operator (JT]) over each of 
these subspaces. Every term in Hamiltonian ([T|) represents a Heisenberg spin exchange that 
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HIDDEN FERROMAGNET IN ABSENCE OF HUND'S RULE COUPLING 
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(2,1 ; 



(2,2) 



FIG. 4: Shown is the low-energy spectrum for 4x4x2 spin-1/2 moments that experience the 
following set of Heisenberg exchange coupling constants: Jo = 0, j\ = 0, J-j 1 > 0, and jj = 
0.3 Ji = J2 ■ The lowest-energy spin-1 state at zero momentum is used as the reference for the 
linear spin- wave approximation. Hereafter, we set h — > 1. 

permutes the Bloch-wave type states. These permutations are stored in memory. Hamilto- 
nian matrix elements that involve Bloch waves composed of configurations of spin up and 
spin down that display absolutely no non-trivial translation invariance up to a spin-flip are 
also stored in memory. They are calculated otherwise. Because the vast majority of Bloch 
waves lie in the first category, the Hamiltonian operator is stored in memory for all practical 
purposes. This speeds up its application on a given state tremendously. The application of 
the Hamiltonian if on a given state is accelerated further by enabling shared-memory paral- 
lel computation through OpenMP directives. Last, we apply the iterative Lanczos technique 
numerically over each subspace labeled by momentum and by spin-flip quantum numbers.— 
The ARPACK subroutine library is employed for this purpose.— 

Figures HI O and [H] show how the low-energy spectrum of the J0-J1-J2 model ([T]) evolves 
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TRANSITION FROM COMMENSURATE SDW INTO HIDDEN FERROMAGNET 
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FIG. 5: Shown is the low-energy spectrum for 4x4x2 spin-1/2 moments that experience the 
following set of Heisenberg exchange coupling constants: jf = 0, Jj 1 > 0, and jj = 0.3 Jj 1 = J^~. 
The Hund's rule exchange coupling is set to Jo = — 1.35 . 

with the strength of the Hund's rule coupling in the presence of off-diagonal frustration. 
The following Heisenberg exchange coupling constants were chosen: j| = 0, Jj 1 > 0, and 
jj = 0.3 Jj 1 = J^- A hidden ferromagnetic state is then expected in the absence of Hund's 
rule coupling, while a cSDW is expected when Hund's rule in obeyed. In the former case 
displayed by Fig. HI notice the coincidence between the linear spin-wave approximation 
about hidden-order, Eq. (|T2|) and Fig. [2j and the present exact-diagonalization results. It 
indicates that the classical prediction of long-range hidden magnetic order holds true for 
spin-1/2 moments per iron orbital. Figure E] displays the spectrum at strong Hund's rule 
coupling, —J = 15 J] 1 , on the other hand. The bulk of the groundstate energy is due to the 
Hund's rule coupling, which is given by 16 Jo = —240 Jj 1 when Hund's rule is obeyed. Notice 
the low-energy spin-1 state at the cSDW wave number (n/a, 0), as well as the softening of 
the spin-1 excited state at the wave number (7c/a,ii/a) that characterizes the Neel state. 
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COMMENSURATE SDW AT STRONG HUND'S RULE COUPLING 
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FIG. 6: The low-energy spectrum of the cSDW is obtained from exact diagonalization of 4 x 4 x 2 
spin-1/2 moments. The following set of Heisenberg exchange coupling constants are used: jf = 0, 
J± > 0, 4 = °- 3 J i = J 2, and Jo = -15 Ji 1 



The dispersion of these spin-1 excited states resembles the spin-wave spectrum predicted 
for the cSDW in the conventional J\-Ji model, where Hund's rule is obeyed.— 1 ^ Figure 
[6] also displays that the groundstate is doubly degenerate, which again is consistent with 
cSDW order over the square lattice. Last, Fig. [5] shows the low-energy spectrum for Hund's 
rule coupling Jo = — 1.35 Jf 1 near a possible quantum phase transition (cf. Fig. [7]). The 
lowest-energy spin-1 excitations at momentum (0,0) and (7r/a, 0) are degenerate, and their 
dispersion resembles the prediction from linear spin-wave theory that is shown by Fig. |2J 
Notably, the lowest-energy spin-1 excitation does not "dip down" at momentum (ir/a, 7r/a). 
Tabled predicts a transition into hidden order at Jo = —0.8 J-j 1 in the classical limit at large 
spin s, however, which suggests that quantum effects renormalize up the critical Hund's rule 
coupling in the physically relevant case of spin-1/2 iron 3d orbitals. 

Figure [7J shows the evolution of relevant magnetic order parameters with Hund's rule 
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It is normalized to its value in the true ferromagnetic state (see text). 



coupling under the previous set of Heisenberg exchange coupling constants. Plotted on the 
vertical axis is the autocorrelation (0(k)± ■ 0(— k)±) of the order parameter 



0(k) ± = ^e ik ^[S J (a)±S i (6)] 



(17) 



over the groundstate normalized to its value in the true ferromagnetic state, (O(0) + ■ 
O(0) + )o = 16 • 17 h 2 . The relation to the ordered magnetic moment per iron atom, /x(±, k), 
given there is obtained from the identity |16^(±,k)| 2 = (2[i B /h) 2 (0(k)± ■ 0(— k)±) that 
is valid over the 4 by 4 square lattice of spin-1 iron atoms. The top graph in Fig. [7| displays 
how the square of the ordered moment for true cSDW order (+) decays with decreasing 
Hund's rule coupling. It notably accounts for the low ordered moment that is observed by 
elastic neutron diffraction in an undoped parent compound to the recently discovered ferro- 
pnictide high-T c superconductors.- Also shown there is how hidden ferromagnetic order is 
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established once cSDW order fades away. Figure [7J then provides evidence for a quantum 
phase transition that seperates a hidden ferromagnetic state at weak Hund's rule coupling 
from a cSDW state at strong Hund's rule coupling. Notice that the putative quantum phase 
transition at J — —1-35 (see Fig. E]) occurs at a Hund's rule exchange coupling that is 
almost a factor-of-2 times larger than the classical prediction J 0c = — 0.8 J-j 1 listed in Table 
HI Last, the bottom graph in Fig. [7] shows the evolution of the relevant order parameters 
in the case of true Neel order with intervening hidden Neel order. Here, we have chosen 
Heisenberg exchange coupling constants jf > 0, J-j 1 = 0, and jj = 0.25 jf = J^, which 
lies at maximum frustration when Hund's rule is obeyed. The common inflection point at 
Jo = — 1.2 j| is consistent with a quantum phase transition at a value of Hund's rule coupling 
that is noticably larger than the classical prediction of Jq c = — jf listed in Table HI 



IV. DISCUSSION AND CONCLUSIONS 



Recent inelastic neutron scattering measurements on the parent compound to ferro- 
pnictide superconductors CaFe2As2 have uncovered well-defined spin-wave excitations in 
the cSDW phase throughout virtually the entire Brillouin zone.^ Instead of softening at 
the wavenumber (jr/a,ir/a) that corresponds to Neel order as predicted by the conventional 
J\-Ji model over the square lattice, they show a local maximum there. Figure [8] displays 
a fit of the spin-wave spectrum obtained by inelastic neutron scattering on CaFe2As2 to 
the prediction for the linear spin- wave dispersion of a square lattice of spin-1 iron atoms, 
(0 + fi_) 1//2 , at the quantum critical point that seperates the cSDW phase from hidden fer- 
romagnetic (Neel) order [see Eq. (fl2"|) and Fig. [2]. The spin per orbital is set to s = 1/2 
and the following Heisenberg exchange coupling constants are used: Jq c = —57.0 meV, 
jj'^ = o, = 115.8 meV, and jj = 43.7meV = These imply a ratio J2/J1 = 0.75 
between the next-nearest-neighbor and the nearest-neighbor Heisenberg exchange constants 
deep inside the cSDW, where Hund's rule is obeyed. Notice that the measured spin-wave 
spectrum terminates before it reaches the predicted Goldstone modes at zero momentum 
that correspond to hidden magnetic order. This is consistent with Fig. |3l which displays 
how the spectral weight of low-energy spin-waves about zero momentum is dwarfed by that 
( 1T5|) of low-energy spinwaves about momentum (jr/a,0) at the quantum critical point that 
seperates the cSDW phase from hidden magnetic order. 
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FIG. 8: The spin excitation spectra of two pnictide materials are fit to the linear spinwave spectrum 
expected at a quantum critical point (QCP) that seperates a cSDW from ferromagnetic (Neel) 
hidden order. [See Eq. (|12p and Fig. [2j] The wavenumber in the direction perpendicular to 
the planes of iron atoms corresponds to antiferromagnetic order in both sets of data. (See refs^ 
ancr^.) 

More recent inelastic neutron scattering measurements on the ferro-pnictide supercon- 
ductor BaFe2_ a: Co 3 ;As2 have uncovered low-energy spin-wave excitations that disperse up 
in energy in an anisotropic fashion away from the wavenumbers that correspond to cSDW 
order;™ e.g. (vr/a, 0). Paradoxically, no evidence for long-range magnetic order exists in 
this ferro-pnictide superconductor. Figure |8] displays a fit of the spin-wave spectrum ob- 
tained by inelastic neutron scattering on BaFe2-a;Co x As2 to the prediction for the linear 
spin-wave dispersion of a square lattice of spin-1 iron atoms, (fi+fi-) 1 / 2 , at the quantum 
critical point that seperates the cSDW phase from hidden ferromagnetic (Neel) order [see 
Eq. ( |T2|) and Fig. |2]. The spin per orbital is again set to s = 1/2, and the following Heisen- 



berg exchange coupling constants are used: Jq c = —105.4 meV, Jj 



IIW 



0, Ji W = 147.2 



18 



meV, and J 2 = 47.3 meV = J^- This fit implies a ratio J2/J1 = 0.64 between the next- 
nearest-neighbor and the nearest-neighbor Heisenberg exchange constants deep inside the 
cSDW. Although the conventional J1-J2 model over the square lattice can also account for 
the linear and anisotropic dispersion of the spinwaves that are observed in superconducting^ 
BaFe2- a; Co :2: As2, it cannot account for the absence of long-range magnetic order. Both the 
semi-classical analysis about hidden magnetic order and the exact diagonalization studies 
on a 4 x 4 lattice described in the previous section indicate, on the other hand, that the 
moment associated with true cSDW order is extremely small at the quantum critical point 
into hidden magnetic order. (See Figs. [Hand El) 

In conclusion, a linear spin-wave analysis of the J0-J1-J2 model Hamiltonian (CQ) that 
describes a square lattice of frustrated spin-1 iron atoms yields a quantum critical point 
as a function of Hund's rule exchange coupling (—Jo) that seperates a cSDW from hidden 
ferromagnetic or hidden Neel order. The observable spectral weight of the Goldstone modes 
associated with hidden order is predicted to be small in general (see Fig. [3]). A strong 
variation between the diagonal and the off-diagonal exchange couplings of the multi-orbital 
Heisenberg model (CQ) across nearest-neighbors alone is sufficient for Hund's rule to give way 
to hidden magnetic order. Recent DFT calculations obtain ferromagnetic direct exchange 
and antiferromagnetic superexchange across nearest-neighbor iron moments in the cSDW 
phase of ferro-pnictide materials.- Their superposition potentially can result in the required 
difference between jf and Jj 1 . Fits of the critical linear spin- wave spectrum to recent 
experimental studies of magnetic excitations in ferro-pnictide materials^ 1 ^ yield a Hund's 
rule exchange coupling constant on the order of 100 meV. Figure [7] displays the dependence of 
the ordered magnetic moments on Hund's rule coupling obtained from exact diagonalization 
of the model Hamiltonian (pQ) over a 4 x 4 lattice, with two spin-1 /2 orbitals per site. 
The top graph implies a critical Hund's rule coupling at the transition between hidden 
ferromagnetic order and cSDW order that is almost twice the classical prediction listed in 
Table [B — J 0c = 0.8 Jj 1 . Correcting for this discrepancy brings the value of the Hund's rule 
coupling closer to the estimate of 700 meV that is obtained for ferro-pnictides by dynamical 
mean field theory— (Note that —Jo is twice as big as the Hund's rule exchange coupling 
that appears naturally in multi-orbital Hubbard models.) 

The success of the theoretical fits to the inelastic neutron scattering data displayed by 
Fig. [8] indicates that a local-moment picture for magnetism in ferro-pnictide materials is 
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valid. Whether magnetic moments in ferro-pnictide materials are itinerant or localized is 
vigorously debated, however. In particular, the metallic temperature dependence shown by 
the electrical resistance in the cSDW phase of ferro-pnictide materials in addition to the 
absence of satellite peaks in inelastic X-ray scattering spectra is consistent with itinerant 
magnetism.— The multi-orbital Heisenberg model Hamiltonian ([[J studied here is also 
potentially consistent with the semi-metallic nature of ferro-pnictide materials, on the other 
hand. Figure [1] indicates that the classical hidden ferromagnetic state remains intact in the 
presence of mobile holes that are prohibited from hopping between orbitals, for example. 
Likewise, the classical hidden Neel state remains intact in the presence of mobile holes that 
are restricted to hop between different orbitals at nearest-neighbor sites of the square lattice 
of iron atoms. Such limiting cases can result in low-energy electronic hole excitations near 
zero momentum, which are consistent with the electronic structure revealed by angle-resolved 
photoemission spectroscopy (ARPES) in the cSDW phase of ferro-pnictide materials.— Res- 
onant scattering of such hole excitations with the quantum-critical spin-waves studied here 
(Figs. [2] and [3]) could also result in low-energy single-particle excitations at momenta near 
those associated with cSDW order. Recent ARPES studies of BaFe2As2 in the cSDW phase 
reveal the existence of low-energy electronic excitations near such momenta.— The mea- 
sured electronic spectrum there forms a Dirac cone, with a Fermi velocity of about 330 
meVA. It is an order of magnitude smaller than in graphene, and it is hence consistent 
with strong electronic correlations.— It also remarkably lies inside the range of spin-wave 
velocities extracted from low-energy magnetic excitations near cSDW wave numbers in the 
doped superconducting compound BaFe2-zCo x As2: 15 v% = 230 meVA and vi = 580 meVA 
(see Fig. [8]). The above observations therefore indicate that injecting mobile holes into the 
states studied here that display hidden magnetic order (Fig. [I]) potentially can account for 
the nature of low-energy single-particle excitations in ferro-pnictide materials. They also 
indicate that a description of the superconducting state of ferro-pnictide materials in terms 
of a doped Mott insulator 7 is possible. 
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